NYU-TH-06/12/10 



Selftuned Massive Spin-2 



Claudia de Rham 1 and Gregory Gabadadze 2 

1 Department de Physique Theorique, Universite de Geneve, 
24 Quai E. Ansermet, CH-1211, Geneve, Switzerland 
2 Center for Cosmology and Particle Physics, Department of Physics, 
New York University, New York, NY, 10003, USA 



Abstract 

We calculate the cubic order terms in a covariant theory that gives a 
nonlinear completion of the Fierz-Pauli massive spin-2 action. The resulting 
terms have specially tuned coefficients guarantying the absence of a ghost 
at this order in the decoupling limit. We show in this limit that: (1) The 
quadratic theory propagates helicity-2, 1, and helicity-0 states of massive spin- 
2. (2) The cubic terms with six derivatives - which would give ghosts on local 
backgrounds - cancel out automatically. (3) There is a four-derivative cubic 
term for the helicity-0 field, that has been known to be ghost-free on any 
local background. (4) There are four-derivative cubic terms that mix two 
helicity-0 fields with one helicity-2, or two helicity-1 fields with one helicity- 
0; none of them give ghosts on local backgrounds. (5) In the absence of 
external sources, all the cubic mixing terms can be removed by nonlinear 
redefinitions of the helicity-2 and helicity-1 fields. Notably, the helicity-2 
redefinition generates the quartic Galileon term. These findings hint to an 
underlying nonlinearly realized symmetry, that should be responsible for what 
appears as the accidental cancellation of the ghost. 



1. Introduction and summary 

An effective field theory for massive spin-2 is motivated by spin-2 QCD resonances 
(glueballs, quark- ant iquark mesons, or their mixtures) which become long-lived in 
the limit of a large number of colors. It is also motivated by massive gravity, and a 
possibility of having dark energy made of "condensate" of massive gravitons. In what 
follows we will be discussing a classical theory of massive gravity on asymptotically 
flat space-time. 

The Fierz-Pauli (FP) Lagrangian [1] is a linear theory describing a massive 
spin-2 state on flat space, without any ghosts or tachyons [2]. Recently, a non- 
linear completion to the FP theory to all orders was proposed in Refs. [3, 4]. 
The gravitational field in this approach is described by an extended metric ten- 
sor g^(x,u), with /j,, v — 0,1,2,3, which is labeled by a continuous dimensionless 
parameter u. The matter fields do not depend on u, but couple to the metric tensor 
g^ u (x) = g^ u [[x,u = 0). The Lagrangian density for the gravitational field alone 
reads as follows: 

C = M&y/gR - £ duy^g - k 2 ) , (1) 

where R is the Ricci scalar of the metric g^ u (x), while k^ u = ^d u g^ u , k = g^k^; all 
indices in the Einstein-Hilbert term in (1) are raised by g^, while those in the second 
term by g^ v . The Z 2 symmetry is imposed on the fields, g^x, u) = g^ v (x, —u). The 
"■u-dimension" is not dynamical since fields have no ordinary derivative terms for 
u > 0, and there are no g uu or g Ufl components of the metric to vary. 

Massive gravity is obtained by requiring that the space-time geometry at u — 1 
is flat. One way to impose this boundary condition is to require that g^x, u = 1) = 
7]^ v = diag(— 1, 1, 1, 1) [3, 4]. Then, the auxiliary dimension can be "integrated out". 
This produces quadratic and nonlinear terms in the resulting effective Lagrangian 
with specially tuned coefficients. In the quadratic order one gets the FP term [3, 4]. 
We will show below that the cubic order result reads as follows: 

C = Ml^gR - **k* (hl„ -h 2 - h^KK + \hhl, - , (2) 

where = g^ u — 77^; all the indices in the mass term in (2) are raised by rf v , and 
the y/gR term in (2) should also be expanded up to the cubic order. 

For small perturbations the quadratic FP term in (2) describes five degrees of 
freedom of massive spin-2. Naively, this would seem enough for perturbative con- 
sistency. Regretfully, the issue is more subtle: to see a potential problem one can 
consider a weak field of a locally-nontrivial asymptotically-flat solution - say a weak 
field of a static lump of matter (below we will refer such backgrounds as local). 
Then, the expression for can be decomposed as 

hfiv = h + 5^, (3) 
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where hr v denotes the weak field of the local background and 5^ denotes small 
fluctuations about it, 1 ^> ^> 5^ u . Substituting (3) into (2) one gets new 
quadratic terms for the fluctuations, h d 55. Since h cl <C 1 these new terms are 
smaller than the 56 terms. Nevertheless, they may destroy the delicate balance of 
the FP term and introduce a ghost. In a generic nonlinear completion of the FP 
massive gravity such a ghost does appear as a sixth degree of freedom [5] . 

In order for the sixth mode not to emerge a very special tuning of the coefficients 
of the quadratic and cubic terms is needed, as was observed in [6] and worked out 
in [7]. What is interesting, is that the theory (1) automatically produces such tuned 
coefficients up to the cubic order! 

There could have been four arbitrary coefficients in front of the five quadratic 
and cubic order terms in (2). In the initial Lagrangian (1) we tuned the coefficient 
between the two terms under the w-integral. As a result, we have automatically 
generated four "good" coefficients in (2) 1 . This hints towards a hidden symmetry 
of (1) which may be responsible for such arrangements. 

One goal of the present paper is to derive the terms presented in (2) from the 
Lagrangian (1). Furthermore, after observing the cancellation of the "ghost terms", 
we will study the remaining nonlinear interactions of the helicity ±2, ±1 and helicity- 
states. 

Thus, we will be looking at the cubic order Lagrangian (2) in the decoupling 
limit where the nonlinear dynamics of all the helicities can be made manifest. This 
method - first used for massive non-Abelian gauge fields in Ref . [8] , and developed 
for massive gravity in [6] - proved to be successful in identifying the presence or 
absence of the sixth degree of freedom in nonlinear theories [6, 9, 7, 10]. The limit 
we consider reads as follows: 

m ->■ 0, Mpi ->■ oo, A 3 = (m 2 Mpi) 1/3 is fixed . (4) 

By taking this limit in (2) we show that: (I) The quadratic theory propagates the 
helicity ±2 modes described by the tensor field h^, the helicity ±1 modes described 
by the vector field A^, and the helicity-0 mode described by the scalar field n - 
all with canonical kinetic terms. (II) The cubic terms with six derivatives, (ddn) 3 , 
which would give rise to the sixth degree of freedom on local backgrounds, cancel 
out automatically due to the special values of the coefficients of the quadratic and 
cubic terms in (2). (Ill) There is a four-derivative cubic term for the helicity-0 mode, 
Dn(dn) 2 , which was first found in the context of the DGP model [11] in Ref. [12]. 
This term is ghost-free for local backgrounds; it is also invariant under the "galilean" 
transformation in the n space, d^n — > d^n+v^, with being a constant four- vector. 
(IV) There are four- derivative cubic terms that mix two helicity-0 fields with one 
helicity-2, such as h((DTr) 2 — (ddn) 2 ) , or two helicity- 1 fields with one helicity-0, such 

1 In general, there exists a one-parameter family of cubic order terms for which the sixth deriva- 
tive terms for helicity-0 do cancel out, (h\ v - h 2 + Cl h u ^h^ - ^^-hh 2 ^ + ^^-h 3 ). The theory 
(1) with the boundary conditions used here gives c\ = —1. From this point of view, (1) had to 
generate only three "good" coefficients up to this order. 
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as 8811(8 A) 2 . All of them are "galilean" invariant, and none of these terms gives 
rise to ghosts on the local backgrounds. (V) If the external sources are ignored (or 
outside of localized sources) all the cubic mixings between the helicity-2 and helicity- 
0, and helicity-1 and helicity-0, can be removed by a nonlinear redefinition of the 
helicity-2 and helicity-1 fields, respectively. What remains is the term, 07r(<97r) 2 . 
Interestingly, the above redefinition of generates also the quartic Galileon term, 
(<97r) 2 ((D7r) 2 — (ddn) 2 [13]. The latter is known to be ghost-free [13] (for more general 
studies of the Galileon, see, [14]). 

That the obtained terms contain the four-derivative cubic term - the cubic 
Galileon \37t(dn) 2 - identical to the one found in Ref. [12] in the context of the 
DGP model, is a hint: The Galileon terms emerge in theories with spontaneously 
broken symmetries [15]. Hence, our findings suggest that there should be an under- 
lying nonlinearly realized symmetry of (1), which is responsible for the cancellations 
of the "ghost terms" . It is conceivable that this symmetry is related to a sponta- 
neously broken 5D reparametrization invariance hiddenly present in the model (1). 
This symmetry should also be helpful in addressing the issue of stability of (1) with 
respect to quantum gravity loops, which has been left open for now. 

So far we have focused on the cubic order. How about an arbitrary n'th order 
terms, (h d ) n ~ 2 55? Making sure that the sixth degree of freedom does not show up 
order- by-order would be a tedious program 2 . However, we are in a better position 
here as the Lagrangian (1) sums up all the polynomial terms in h^ v . Then, the 
presence of a ghost could be seen by calculating the exact Hamiltonian. Boulware 
and Deser (BD) [5] have shown that the sixth degree of freedom would in general 
appear in nonlinear massive gravity due to the loss of the Hamiltonian constraint: 
This leads to a Hamiltonian term that is proportional to positive powers of the 
canonical momenta, but is sign indefinite, hence, signaling the presence of a ghost. 

Therefore, the absence of the BD term would be a good indicator that the sixth 
degree of freedom is not present 3 . The Hamiltonian for (1) was calculated in [3], 
where it was shown that the BD term cancels out. Hence, one should expect that 
the BD ghost does not appear in the order- by-order expansion of (1). Moreover, in 
Ref. [4] the decoupling limit of the theory was considered and it was shown to all 
orders that the leading terms arising at the scale A < A 3 , that could give rise to 
ghosts, cancel out 4 . 

There have been proposals in the literature to obtain the theory of massive spin- 
2 via a dynamical condensation mechanism (for recent works see, e.g., [17, 18] and 

2 Since in each consecutive order these terms have smaller and smaller coefficients, one could 
try to show that after a certain order the ghost appears only above a certain UV cutoff of the low 
energy theory. 

3 It is not a guarantee, however, of the absence of other types of ghosts that may be present in 
a theory for some other reasons (e.g., introduced by hand). Here, we're focusing on the ghost that 
may appear in massive gravity due to the nonlinear interactions. 

4 Ref. [7] concluded that no linear combination of the quartic order terms in h^ v can give a 
theory for n that would be ghost-free. This issue will be revisited in our forthcoming paper [16], 
with a different conclusion. Until then, we ignore the explicit quartic and higher terms in h^ v . 
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references therein). It would be interesting to see whether the cubic terms in these 
models can also automatically give rise to ghost-free structures for the n field. 



2. Integrating out the auxiliary dimension 

The goal here is to calculate order-by-order the w-dependence of the extended metric 
gn U (x,u), then substitute it back into (1), and integrate the latter w.r.t. u. This 
should give the effective Lagrangian written in terms of only. To fulfill this goal 
we introduce the notations 

~ 9lxv (x, u) = Vfll/ + if <£> (x,u) + H$ (x,u) + H$ (x, «) + ... , (5) 

where H$(x,u), H^J(x,u), ... are perturbations in the corresponding order. Since 
g^ v (x,u = 0) = g^„{x) = r]^ + h^x), we require that 

H${x,u = 0) = h^x), H$(x,u = 0) = H$(x,u = 0) = 0. (6) 

Furthermore, following Refs. [3, 4] we impose the boundary condition on the ex- 
tended metric g^x^u = 1) = r]^ u , which guarantees that the theory reduces to FP 
massive gravity in the quadratic approximation. This boundary condition implies 
that HjiJ(x,u = 1) = H$(x,u = 1) = H$(x,u = 1) = 0. Having set these, the 
calculation of the li-dependence of g^ u (x,u) is well defined. 

The expression for the mass term in (2) in terms of H$(x, u), H^J(x, u), up 
to cubic order, is straightforward to obtain: 

_M|^ 2 £i du ^ duH ^ )2 + 2duH (i JduH m _ 2HM du H%d u HM - {d u H^ f 

-2d u H^d u H^ + 2d u H^H$d u H$ + \^\d u H^f - l -H^\d u H^f), (7) 

where we used the fact that the metric is a Z 2 symmetric function of u, and set the 
integration limits from to 1. Also, in (7) and in what follows, we use simplified 
notations with all lower-case indices contracted by the flat space-time metric. 

To find the w-dependence of H$(x,u), Hf$(x,u), ... we integrate the equations 
of motion. For this we vary the action (1) w.r.t. g^ u (x,u). The resulting equations 
for u = + and < u < 1 read respectively as follows: 

Gfj, v - m 2 (kpv - g^k) = T M ,/(2M| 1 ) , (8) 

and 

d u [y/i {kr v - kn] = (^ 2 - + ^ - «) . (9) 

The latter is the equation that determines the "w-evolution" of the extended metric. 
With the boundary conditions specified above it is straightforward to solve (9) for 
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H^J(x,u), HlJ(x,u), .... It turns out that only the solution for Hj^J(x,u) is needed 
at the cubic order. The latter reads: 

ff«(s,u) = (l-u)V(*)- (10) 

This expression for H$ (x, u) guarantees that the two terms in (7), which contain the 
function H^J(x,u), integrate to zero due to the boundary conditions. Substituting 
the expression (10) into (1) and integrating it w.r.t. u we obtain 

-^p- fa -h*- K kk + \hhl, - \h^j , (ii) 

which is the mass term presented earlier in (2). 



3. Extracting the longitudinal mode 

We begin by rewriting the Lagrangian density (2) in a manifestly covariant form. 
For this, following Ref. [6], we introduce a covariant tensor field H^ v , which is 
related to g^ u as follows: 

dr{x)d4>\x) 

H IHf = g lw -r hb ^ Qxv , (12) 

where <p a (x) are just four scalars; a, b = 0, 1, 2, 3, and r] a b = diag(— 1, 1, 1, 1) is the 
flat metric on the field space of the scalars. This construction guarantees that H^ v 
transforms as a covariant symmetric rank-2 tensor. 

Furthermore, it is convenient to decompose the scalars as 4> a (x) = x a — ir a (x), 
where x a = S^x 11 . This decomposition specifies that under the general coordinate 
transformations, x M — > x^ + (^{x), the fields ir a (x) transform as 7r a — > n a + S^C 1 *- 
Using the above definitions, we can easily find the expression for the tensor in 
terms of = h^/M^V- 

„ V , + dyVjx _ d^V a d v V a 

where = 5^Tr a Al, is a field that shifts as — > + rj^^Al under the general 
coordinate transformations (the index contraction in (13) is done with T]^). 

Then, the Lagrangian density (2) can be written in terms of the covariant tensors 
only. It reads as follows: 

C = M^yfgR - (h% -H> + H;h:H» - \hH% + , (14) 

where all the indices are raised with the metric tensor g^ u (as a consequence, the 
signs in front of the cubic terms in (14) flip as compared with (11)). In the cubic 
order (14) reduces to (2) after using the gauge fixing condition = 0. 
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Moreover, the external sources/fields couple to g^ u : For instance the Lagrangian 
for a scalar field tp coupled to gravity would read as follows: 

= Ivgi-g^d^ - 2vm. 

We will not write explicitly the couplings of g^ v to the external sources/field below, 
but will keep them in mind (see, discussions below). 

To turn to the decoupling limit we decompose the vector field as follows: 

V ll = mA ll + d lt ir, (15) 

where both and ir are kept finite in the limit (4). The field A^ will end up 
encoding the helicity ±1 states, while the field ir will describe the helicity-0 state. 

Next we calculate the decoupling limit of the Lagrangian (14). This is done by 
substituting (13) and (15) into (14), and taking the limit (4). This procedure - 
valid for fields that decay fast enough at spatial infinity - requires some care: we 
introduce an infrared regulator of the theory, say a large sphere of radius L >> 1/m, 
and take the radius to infinity, L — > oo, before taking the limit (4). This hierarchy 
of scales enables us to put all the surface terms to zero before taking the decoupling 
limit. 

Once the above procedure is adopted we find the following remarkable properties: 
(1) All the terms containing six derivatives and three helicity-0 fields, such as (<9 2 7r) 3 , 
that come suppressed by the scale A 5 = (Mpim 4 ) 1 / 5 <C A 3 , cancel out up to a total 
derivative [4]. (2) The quadratic terms in A form the Maxwell term, while all the 
terms that are linear in A and quadratic in ir, such as dAd 2 nd 2 Tr, which would 
be suppressed by the scale A4 = (Mpim 3 ) 1 / 4 <C A3, also cancel out up to a total 
derivative. (3) The only terms that survive are those suppressed by the scale A3. 

In this section we focus only on the helicity-2 and helicity-0 modes, while the 
terms with the helicity-1 field will be ignored until the next section, where they are 
shown to be harmless. 

The remaining terms, after the conformal transformation h^ u = + f]^ u ii that 
diagonalizes the quadratic action, read as follows: 

1 h FV*"#h 1 3 ^ 1 3D7r(<yr) 2 

h((d^d u 7i) 2 - (Pit) 2 ) /y (d^d u nDn - d^d a nd u d a Tr) 
4A| 2A| 

Here, all the indices are raised using the flat space metric and we do not distinguish 
between the upper and lower cases. £ denotes the Einstein operator that is related 
to the linearized Einstein tensor G^ v as follows: 

S^h afS = G» v = -i(D/^ - d"d a h a - d u d a h» + - if v Uh + rTd a d p h ap ). 
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All the terms in (16), up to total derivatives, are invariant under the "galilean" 
transformations of the 7r field, d^n — > d^n + t> M , where is some constant four- 
vector. Moreover, none of the cubic terms in (16) produce ghosts on any local 
background. Indeed, the last term in the first line of (16) is identical to the one 
found in the decoupling limit of the DGP model [12] (see, also [19] for related 
discussions). This term is known to give rise to the equations of motion that have 
no more than two derivatives acting on each fields [12]. To see that the rest of the 
terms in (16) are also safe we rewrite them as follows: 

h ((^tt) 2 - (Dtt) 2 ) = h (2d 2 nAn - 2(d d, j n) 2 - (Avr) 2 + (c^-tt) 2 ) , (17) 

and 

h^ u (d^d v nOn - d^d a nd u d a ir) = h 00 (d^nAn - (d djir) 2 ) + 
2h j (d djirAir — d d k irdjd k ir) + hij (didjirOn — did^ndjO^n) . (18) 

These have at most two time derivatives, and will not produce any terms with 
more than two time derivatives in the equations of motion. Importantly, the term 
multiplying hoo in (16) has no time derivatives, only spatial ones, showing that Jiqo 
remains a Lagrange multiplier at the cubic order in the decoupling limit. 

The expression (16) is invariant under the gauge transformations, Sh^ = —d^— 
d u (^, due to the cancellation between the two nonlinear terms in the second line. In 
other words, the Bianchi identity is automatically satisfied for the tensor equation 
that follows from (16) by varying it w.r.t. h^. 

If we ignore coupling to external sources, one can simplify further the Lagrangian 
(16) by the following nonlinear transformation 

V = Ku + ^3 d v nd » n ■ ( 19 ) 
The resulting Lagrangian reads: 

= -h'^h'a, + + + ••■ (20) 

As emphasized before, in the present theory it is the field that couples to external 
sources/fields. In the linearized theory the linearized source/field stress-tensor T^ u 
couples as h^ u T^ u /M P i, which after the conformal transformation reads as follows 
(h 'pJTpv + TtT)/M-p\. These couplings, e.g., for a static source, are held fixed and 
finite in the decoupling limit (i.e. T^/Mpi is fixed to be finite [12]). However, 
the diagonalization of the nonlinear terms performed by (19) would generate an 
additional nonlinear coupling d^d^T^/ (M?\h%). Hence, to avoid complications 
with the additional nonlinear couplings and the modified light-cone, it is better to 
think of the "decoupling" limit Lagrangian (16) in which the helicity-2 and helicity-0 
field mix at the cubic order. For certain sources - such as static ones - the additional 
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coupling d^nd^T^/^MpiAf) is zero, and therefore using the Lagrangian (20) in 
which the helicity-2 and helicity-0 are truly decoupled may be more convenient. 
Interestingly, the field redefinition (19) in (16) generates in (20) the term 

2d f ,nd u 7r(d^d l3 7rd l/ d /3 7T - d^nDn) + (<9 M 7r) 2 ((D7r) 2 - (d^nf) 

16AI • 1 ] 

This is exactly the quartic Galileon introduced in [13] as a ghost-free quartic order 
term giving rise to two derivative equations of motion. Note that the coefficients of 
the cubic Galileon in (20) and the quartic Galileon in (21) are related to each other 
since they both originate in the cubic mixing terms of helicity-2 with helicity-0 in 
(16). In the present paper we ignore the quartic order terms in h^, but in general, 
depending on the coefficients of these terms, the quartic Galileon (21) may or may 
not cancel. These issues will be discussed in detail in [16], where it will be shown 
that in a general quartic-order theory, as soon as the cubic Galileon is present in 
(20), we are also bound to generate either the quartic Galileon, or a quartic mixing, 
or even both together accompanied by the quintic Galileon. 

In practice, the quartic and higher order terms are negligible at large scales 
where the tensor-scalar gravity sets in [20], while they become as important as the 
quadratic and cubic contributions around the Vainshtein scale [21] (see also [22]). 
The mixing terms will be essential to address the issue whether the full theory admits 
superluminal propagation, as the pure Galileon terms do, (see Refs. [23, 13]). 



4. Extracting the helicity-1 modes 

The Lagrangian for the vector field in the decoupling limit is obtained by substituting 
(13) and (15) into (14), and taking the limit (4). 

As mentioned before, the quadratic terms in A form the Maxwell term. Moreover, 
terms linear in A and quadratic in n, (e.g., dAd 2 nd 2 ir), that would be suppressed 
by the scale A 4 = (Mpim 3 ) 1 / 4 <C A 3 , cancel out up to total derivatives. The terms 
that survive in the decoupling limit are suppressed by the scale A 3 . These mix 
two helicity-1 fields with helicity-0, dAdAd 2 ir. The resulting Lagrangian reads as 
follows: 

C A = --F 2 U -^j^(2d^A a d v A a + 2d a A IM d a A v + 8d IM A a d a A v -12d a A a d li A l/ ) 
- ^ {-{daAp) 2 - 5(d„A v d u A„) + 6(d a A a ) 2 ) . (22) 

As before, all the contractions are by r]^ u , and no distinction is made between the 
lower and upper-case indices. 

The above expression is invariant under the internal galilean transformations. 
However, the gauge invariance of (22) is not immediately obvious. We leave it to 
reader's pleasure to show that (22) reduces, up to a total derivative, to the following 
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Lagrangian: 




lp 2 - — 

4 2A| 



- r]^D)(d a nF ai/ ) . 



(23) 



The above expression is invariant, up to a total derivative, w.r.t. the gauge trans- 
formations of the vector field, — > + d^Xi where x is the gauge parameter. 

Furthermore, the nonlinear terms in (23) can be removed by the following non- 
linear field redefinition: 



After this transformation, and up to quartic terms, we are left with the Maxwell 
Lagrangian. Therefore, (22) describes helicity ±1 modes, and does not give rise to 
ghosts at cubic order on any local background. 

5. Brief comments 

Since the summary of our main results has already been given in Section 1, we end 
this work with a few technical comments. 

(i) In Section 3 we calculate the decoupling limit of (14) using the method of 
Ref. [6]. The method of Ref. [10], although similar to that of [6], differs from 
it slightly and follows more closely the Stuckelberg method for the gauge fields. 
Furthermore, we have checked that in the approach of Ref. [10] all the six-derivative 
cubic terms cancel out, and four derivative ones remain. The mixing terms for 
helicity-2 and helicity-0 are present, and take a somewhat different form, but satisfy 
gauge invariance and yield automatically the Bianchi identities. One can also show 
that those mixing terms are reducible by a nonlinear field redefinition to the ones we 
obtained here, and are also removable by yet another field redefinition (if external 
fields/sources are ignored) at the expense of generating the higher order terms. 

(ii) It is straightforward to see that the special coefficients in (2) also play a 
role in the construction of the Hamiltonian. Indeed, introducing the standard ADM 
decomposition [24] with the lapse N and shift Nj we find that the cubic order 
Hamiltonian of (2) is linear in 5N = N — 1. Hence the Hamiltonian constraint is 
maintained in this order 5 . On the other hand, Nj enters quadratically and is alge- 
braically determined, as required for massive spin-2. Note that the decoupling limit 
considerations do not yet guarantee positive semidefiniteness of the Hamiltonian of 
the full theory, and this has to be addressed separately (see discussions in [3]). 

(iii) Ref. [7] showed that 5N does get quadratic terms in the quartic order of a 
general massive theory (see also [5, 25]). This dependence, however, may come in 
a special combination with Nj that allows to preserve the Hamiltonian constraint 
and avoid the BD term, see more on this in Ref. [16]. 

5 Note that hoo enters quadratically, but this does not prevent the Hamiltonian to be linear in 
SN, which is the right variable at the nonlinear level and away from the decoupling limit. 
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(d a 7iF afJ: ) . 



(24) 



" 2AI 
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